Scanning tunneling spectroscopic evidence of crossover transition 
in the two-impurity Kondo problem 
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We calculate the differential conductance (dl/dV) corresponding to scanning tunneling spec- 
troscopy (STS) measurements for two magnetic atoms adsorbed on a metal surface with the aid 
of the numerical renormalization group (NRG) technique. We find that the peak structure of the 
dl/dV spectra near the Fermi level changes gradually as a function of the adatom separation and 
the coupling between the adatoms and the metal surface conduction band. When the coupling 
becomes small, the peak disappears and, instead, a dip structure appears near the Fermi level. This 
dip structure is the manifestation of the strong antiferromagnetic correlation between the localized 
spins. The gradual change of the dl/dV structure from a peak structure to a dip structure originates 
from the crossover transition in the two impurity Kondo problem. 

PACS numbers: 73.20.-r, 75.30.Hx, 75.75.+a, 07.79.Cz, 75.40.Mg 



I. INTRODUCTION 

The two-impurity Kondo problem has long been ex- 
tensively investigated. Generally, the low temperature 
physics depends on the ratio between the RKKY coupling 
constant Jrkky and the one-impurity Kondo tempera- 
ture Tk In the limit of strong ferromagnetic RKKY 
coupling (—Jrkky >> Tk), a two stage Kondo effect 
occurs: First the localized spins are partially compen- 
sated as temperature T decreases, and then completely 
suppressed as T goes to zero. In the limit of strong an- 
tiferromagnetic RKKY coupling (Jrkky >> Tk), the 
localized spins are locked into a spin-singlet state (the 
antiferromagnetic region), and the Kondo effect plays a 
minor role. In the Tk >> | Jrkky | region (the Kondo 
region), the RKKY interaction plays a minor role. 
The scenarios in the corresponding limits are very rea- 
sonable, but the physics in the region at intermediate 
values of Jrkky is nontrivial. Studies carrying out nu- 
merical renormalization group (NRG) diagonalization on 
the electron-hole (e-h) symmetric Hamiltonian found a 
critical point which separates the Kondo region and the 
antiferromagnetic region However, it is now appar- 

ent that the critical point results from the e-h symmetry 
of the model Hamiltonian0,Q. In general, e-h symmetry 
is broken (e.g., energy dependence in tunneling matrix el- 
ements) and the quantum phase transition is replaced by 
a crossover transition Q . 

The magnetic atom dimer on a nonmagnetic metal sur- 
face is a classical example of a two-impurity Kondo sys- 
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tern. In this system, the strength and sign of the RKKY 
interaction can be adjusted by changing the adatom 
separation. Using STS, the Kondo effect can be ob- 
served through the sharp peak structure near Fermi 
level, which corresponds to the Yosida-Kondo singlet 
0, i, ©, M, El E El 0, El • Thus, we would expect to 
observe the interference between the Kondo effect and the 
spin ordering effect of the RKKY interaction through the 
STS spectra. Experimentally, the STS spectra of a Co 
dimer on a Cu(100) surface vary with the change of the 
adatom separation [171] ■ which is expected to result from 
the two impurity Kondo problem. However, theoretical 
studies specific to STS measurements are rare [u3. [THI . [l6| . 
In our earlier studies, we confirmed that the ferromag- 
netic (antiferromagnetic) interaction tends to sharpen 
(broaden) the dl/dV peak structure [13, fig . On the 
other hand, we find from dl/dV calculations based on e-h 
symmetric Hamiltonian that such a broad peak does not 
exist. These results indicate that the observed broad STS 
spectra are related to the crossover between the Kondo 
region and the antiferromagnetic region. 
Both Tk and Jrkky strongly depend on the ratio be- 
tween the Coulomb interaction U and the coupling be- 
tween the adatom and the metal surface conduction 
band, V. As schematically shown in Fig[TJ Tk decays 
much faster than Jrkky as T decreases. With small F, 
we would expect that the antiferromagnetic RKKY in- 
teraction becomes dominant. In this situation, adjusting 
the RKKY interaction by changing the adatom separa- 
tion, we can trace the crossover between the Kondo re- 
gion and the antiferromagnetic region. For this reason, 
we calculate the dl/dV spectra for several values of T and 
the adatom separation to show how the crossover transi- 
tion in the two impurity Kondo problem can be observed 
through the STS spectra. 
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FIG. 1: Diagram of the smooth crossover in the two impu- 
rity Kondo effect. T is the temperature. Tk is the Kondo 
temperature and Jrkky is the RKKY coupling constant. 



II. MODEL AND METHOD 



As experiments show[Tl|, [l8j], by covering the metal 
surface with a decoupling layer, T was suppressed, which 
leads to a decrease in Tk- In the present study, we con- 
sider the model system shown in Fig(2J 
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FIG. 2: STS observation of a magnetic dimer on a nonmag- 
netic metal surface. Distances between the adatoms (a) , the 
tip apex and the metal surface (z p ) are given in An gstroms 
[A]. STM tip is placed directly above atom 1. 
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Here d icr , c ka and c prT correspond to creation operators for 
adatom d electrons, metal surface conduction electrons, 
and tip electrons with spin a, respectively, = d\ a di n . 
Ed, Ek and E p correspond to the kinetic energies of 
adatom d electrons, conduction electrons and tip elec- 
trons, respectively. Adatom index i = 1,2. k corresponds 
to the metal surface electron wavenumber, and p corre- 
sponds to eigenstate quantum number of the tip elec- 
trons. T p di, W p k, and Vkdi correspond to the tip-adatom, 
tip-surface, and adatom-surface electron tunneling ma- 
trix elements, respectively. U gives the on-site Coulomb 
repulsion on adatoms. We approximate the coefficients 
in the Hamiltonian fTJ) as 



Vkdi 
W k p 



' pdi 



Vbexp(zk • a/2), V kd2 = V exp(-ik • a/2), (2) 
W exp(-(z p - R s + i? )/A)exp(-z(k • r p )), 

(3) 



ipdi{r P - Rj) 

J 0-7 — 7R — T^TV — -LQ^di, 



^i(Ri + R) 



and 



-0 di (r p - Ri) 



(4) 



(5) 



^(Ri + R) 

ipdi corresponds to the d electron orbital of adatom i. 
We approximate the tip apex as a nonmagnetic metal 
sphere whose center is positioned at r p = (a/2, 0, z p ) 
with radius Rs. R = (0,0, Rs) and a = (a, 0,0). Wo 
and To correspond to the values of the tunneling matrix 
elements when the tip is in contact with the surface or 
adatoms, respectively. Vq is the tunneling matrix clement 
between a localized d electron and a metal surface in the 
single impurity case. 

Using non-equilibrium Green's function method[9LIT9j|. 
the electron current from the STM tip to the surface can 
be written as 



Magnetic atom 2 is located at a distance of a A with 
respect to magnetic atom 1. We assume that the center 
of atom 1 is located at Ri = (a/2, 0, R m ) A and that of 
atom 2 is located at R 2 = (—a/2,0, Rm) A. R m is the 
radius of adatom. We set R m = 0.4 A @. The tip apex is 
at a distance of z p from the metal surface, directly above 
atom 1. R is the thickness of the decoupling layer and 
we set R a — 5A. The corresponding model Hamiltonian 
is given by 



1 = T^J ^(A-/p){ 2 ^ 2 ft + ^K3GI lff 

+ 4irT^di<5>d2Pp%G r 12a - 2TT 2 p kP pW^e~ 2( - z "~ R ^/ x 

- 47rT ($«ii J ul0 (n) + $ d2 Jvw(r 2 ))Pp^G r lla 

- 4nT ($ d iJ vw (r 2 ) + $d2Jvw(ri))p P $lG r 12u 

- 27r(J^(n) + J 2 w (r 2 ))pp^Gl la 

- 4:TrJ vw (r 1 )J vw (r 2 )pp'3G r 12(J \. (6) 



3 



From Eq. ([6]) we can then obtain the corresponding dif- 
ferential conductance dl/dV. In Eq. ©, the retarded 
Green's function Gy a of the surface system corresponds 
to that of the two-impurity Anderson model, fk and f p 
give the Fermi distribution functions for the surface and 
tip, respectively, pk and p p give the density of states 
of the conduction electrons and the STM tip electrons, 
respectively. We define $> d i and J vw (r) as follows - 

Jvw(r) = TrpkJo(k F r)VoWoexp(-(zp - R s + R a )/X). 

(7) 

A is the decay constant of the surface electron wave func- 
tion. Jo is the Oth order Bessel function. Using these 
parameters, T is defined as T = ttpuVq. At T = K, 
dl/dV can be rewritten in the following simplified form 



dI/dV(eV) oc A n <SG r U(7 {eV) 
+ B n ^G r U(7 (eV) 



A 12 -sG r 12a (eV) 
■ B 12 m r 12a (eV) (8) 



Coefficients such as An, Bn can be derived from the co- 
efficients of the Green's function in Eq. ©, i.e., 

A n = -2irT%{$ dl + $2 2 ) ppj 

+ M-On) + JL(r2))p P , (9) 

A 12 = -4:TrTfa dl <S> d2 p p +iTrJ vw (r 1 )J vw (r 2 )p p ,(10) 

B n = iTTT (^ d iJ vw (n) + ^ d2 J vw (r 2 ))p p , (11) 

B 12 = 4nT (<I> dl J vw (r 2 ) + <f> d2 J vw ( ri ))p p . (12) 

The third and fourth terms are related to the Fano ef- 
fect and makes the dl/dV spectra asymmetric. We show 
several values of each coefficient eq.® in the case of 
r = 0.022eV at several adatom separation in Table U 



TABLE I: Values of coefficients in Eq.® in the case of T = 
0.022 eV at several adatom separations (a). 



As can be seen from TablelH in eq.©, Ai 1 ( ^G r lla {eV) 
is dominant and other parts play only minor roles. This 
means that the decoupling layer suppresses the Fano ef- 
fect. To derive the relevant Green's function, we adopt 
the NRG technique ||, EH [H [H, SI. In the present 
study, we transform H^ 2 in eq.([T]) into two semi- infinite 
chains form so as to be suitable for NRG calculation 



HA2 - ^ e nqlnqfnq + ^ {tnqP nq fn+lq + H.C.) 
nq—i: nq—i: 

ndi-\n d ii 

q=± 1=1,2 



2DT 



Y,(v^riqd q + h.c). 

<r=± 



(13) 



Here, f nq corresponds to the nth site of the conduction 
electron part of the chain (Wilson chain) with the par- 
ity q. q = +, — denotes the even and odd parity states, 
respectively. d + — dl ^ 2 and <2_ = dl ^ 2 ■ D gives the 
conduction electron band width. e nq and t nq are the cor- 
responding matrix elements for the Wilson chain. When 
the metal surface conduction band is that for a two di- 
mensional free electron, 



J-D 

1 ± J (/ea) 



7±(e) = 
7± = 



D 



— D 



7±(e)de. 



(14) 
(15) 
(16) 



a (A) 


An (eV) 


Am (eV) 


Bn (eV) 


B 12 (eV) 


5.0 


-6.6225xl0~ 4 


-1.6792xl0~ 9 


5.3360xl0~ 7 


1.7792xl0~ 7 


6.0 


-6.6225xl0 -4 


6.8125xl0~ 10 


5.3360xl0~ 7 


6.6775xl0" 8 


7.0 


-6.6225xl0~ 4 


2.9035xl0~ n 


5.3360xl0~ 7 


-3.5216xl0~ 8 


8.0 


-6.6225xl0~ 4 


-4.6240xl0" n 


5.3360xl0~ 7 


-1.1929xl0 -7 


9.0 


-6.6225xl0~ 4 


-7.1997X10- 11 


5.3360xl0~ 7 


-1.7886xl0~ 7 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

The two dominant parameters Tk and Jrkky 
strongly depend on the value of T. Tk decays much faster 
than Jrkky as T decreases. The values of T (therefore 
Tk) depend on the thickness, the surface condition and 
the type of the decoupling layer. In the present study, we 
set the thickness of the decoupling layer to 5 A [11| . In 
such system, the value of Tk ranges from 3 to 6K. Thus 
we calculate the dl/dV spectra with T = 0.025,0.022 
and 0.02eV at several adatom separations. (From [2~i| 
corresponding Tk is estimated at 6.64, 3. 07, and 1.63K 
respectively) In Fig|3j we show the calculation results for 
the dl/dV spectra. At a — 5.0A and a = 9.0A, there is 
a sharp peak structure near the Fermi level. With these 
adatom separations, the RKKY interaction is weak anti- 
ferromagnetic and the Kondo effect is dominant [l5j . The 
sharp peak corresponds to the Yosida-Kondo resonance. 
With our settings, the antifcrromagnetic RKKY interac- 
tion becomes largest around a = 7.0A. As the adatom 
separation becomes close to 7.0A, the peak structure 
changes gradually. When T = 0.025eV, the dl/dV spec- 
tra broaden as a becomes close to 7.0A. However, when 
r = 0.022eV, a dip structure appears near the Fermi 
level. This dip structure develops to a deeper one as T 
decreases. 

We find that the "parity splitting" of the single elec- 
tron excitation spectra is the origin of the dip structure 
in the dl/dV spectra. As shown in eq. (fT3"|) . the opera- 
tors in the model Hamiltonian are indexed by the parity 
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FIG. 3: dl/dV calculation results at several adatom sepa- 
rations with a) T = 0.025eV {T K =6.64K), b)T = 0.022eV 
(T K =3.07K), and c) F = 0.02eV (T K =1.63K). The STM tip 
is placed at z v = 4.0A. We set R s = 3.0A and Wo = To = 
0.02eV. 



q. Thus, the obtained physical properties such as sin- 
gle electron excitation spectra are divided with respect 
to the parity. The dl/dV spectra are proportional to 
the average of the single electron excitation spectra of 
adatom electrons on each parity channel (p+,p-). As 
shown in Figdl p + and p_ have different peak positions. 
The asymmetry of the spectra results from the difference 
in the coupling between the adatom and conduction band 
in each parity channel (See eq. ([P5)) and eq. (fl~5)) ). Previ- 
ous studies shows that the parity splitting is one of the 
signatures of the smearing out of the critical point due to 
the breaking of the e-h symmetry [1, 0] . In the present 
study, based on the model Hamiltonian (fT]), the tunnel- 
ing matrix element between the localized d-electrons and 
the metal surface conduction electrons has energy depen- 
dence, which breaks the e-h symmetry. 

In order to investigate the electron state, we 
plotted the flow of low lying many-particle energies 
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FIG. 4: Single electron excitation spectra of adatom electrons 
for each parity channel. 



in the NRG calculation with several adatom separa- 
tions (Fig[5]). The states are labeled by the quan- 
tum numbers total charge Q, total spin S, and to- 
tal parity P. Here we show the flows of the lowest 
state for (Q,2S,P) = (0,0,1),(1,1,1),(-1,1,1),(1,1,-1),(-1,1,- 
1), (0,2,1) and (0,2,-1). The Q=0,2S=0,P=1 state is the 
ground state. Q=±l,2S=land P=±l states are the sin- 
gle electron (hole) excited states. When a = 5.0A, as 
N increases, the energy of (±1, 1, ±1) state decreases 
relatively. This indicates the energy gain from the for- 
mation of a singlet between the adatom localized spin and 
excited state of the conduction electron, i.e., the Yosida- 
Kondo singlet. However, when a comes close to 7.0A, 
the energy of the single particle excited states do not de- 
crease so much but keep considerable value in large N. 
From the calculation result for the spin correlation func- 
tion between adatom localized spin(Fig[6]), we conclude 
that this large energy of single particle excitation states 
originates from the interruption of the Yosida-Kondo sin- 
glet formation by the strong antiferromagnetic correla- 
tion between the localized spins. In our model, the e-h 
symmetry breaking connect the antiferromagnetic region 
and the Kondo region continuously and the ground state 
would be the hybrid of the localized spin-spin singlet 
state and the Yosida-Kondo singlet state. As the anti- 
ferromagnetic correlation between localized spin becomes 
large, the localized spin-spin singlet state would become 
dominant and suppress the excitation at the Fermi level, 
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FIG. 5: Flows of the low- lying many-particle energy levels in 
odd iterations. N + 1 corresponds to the iteration step num- 
ber. In this calculation we calculate at a — 5.0, 6, 0, 7, oAa) 
r = 0.025eV b)r 0.02eV. The states are labeled by the quan- 
tum numbers total charge Q, total spin S, and total parity 
P. 



which results in the dip structure of dl/dV. Thus, the dip 
structure in dl/dV would enable us to detect the strong 
antiferromagnetic correlation between localized spins. 

We can also discuss about the fixed point Hamil- 
tonian and the phase shift estimated from that. As 
shown in FigjSl the energy levels converge to some values 
for large N. This means that the Hamiltonian becomes 
close to the fixed point Hamiltonian of rcnormalization 
transformation [2fjl l2ll [23[. In this calculation, we can 
define the fixed point Hamiltonian as follows; 

N-l 
q=±,n=0 

+ K q flh) (17) 

q = ± is the parity index. The second term with K q 
is related to the influence of potential scattering in each 
channel. The strength of the potential scattering K q and 
the phase shift 8 q have the following relation: 
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FIG. 6: Spin correlation function as a function of adatom 
separation. 
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FIG. 7: Phase shift as a function of adatom separation, (in 
units of 7T radians) 



We estimate K q in each conduction channel by compar- 
ing the many-particle energy of first excited state in the 
last iteration of NRG calculation with that of Eq. (fT7j) . 
The phase shift is related to the number of electrons 
and holes which are virtually bounded by localized spin 
- i.e., the phase shift is a barometer of the existence 
of the Kondo effect. If 8± — tt/2, it means that one 
electron-hole pair is bounded to each localized spin and 
quenches it individually (i.e,. two Yosida-Kondo singlets 
are formed). On the other hand, 6± = indicates that 
a localized spin-spin singlet is formed 0. In FigH we 
show the result of the phase shift calculation as a func- 
tion of adatom separation with several values of T. 
The phase shift changes gradually and has an intermedi- 
ate value between ±7r/2 to (Fig. [7|). The smallest value 
of S± becomes close to as T decreases. This result also 
indicates that the dip structure around a = 7.0A orig- 
inates from the dominance of the spin-spin singlet state 
in the ground state. Though the values of 5± at each 
adatom separation are different by T, the change of S± 
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is smooth in all cases. If a critical point separates the 
Kondo region and the antiferromagnetic region, the pos- 
sible values of 5± are only or ±7r/2 0, [y, [25| . The 
intermediate value of 6± would result from the crossover 
transition and supports the conclusion that there is no 
critical point in the system of a magnetic dimcr on a 
metal surface. 



IV. SUMMARY 

To investigate how the transition between the Kondo 
effect dominant region and the antiferromagnetic RKKY 
interaction dominant region can be observed through 
scanning tunneling spectroscopy (STS), we calculate 
the differential conductance (dl/dV) corresponding to 
STS measurements for two magnetic atoms adsorbed 
on a metal surface with the aid of the numerical 
renormalization group technique. We find that the peak 
structure of the dl/dV spectra changes gradually as a 
function of the adatom separation and the coupling (r) 
between the adatoms and the metal surface conduction 
band. When T becomes small, the peak disappears 
and a dip structure appears near the Fermi level. This 
dip structure originates from the parity splitting of the 
single electron excitation spectra and the manifestation 
of the strong antiferromagnetic correlation between the 
localized spins. The result of the phase shift calculation 
supports the conclusion that there is no critical point 
between the Kondo effect dominant region and the 
antiferromagnetic RKKY interaction dominant region 
but a crossover transition connect these regions. 



In conclusion, we show that the crossover transition from 
the Kondo region to the antiferromagnetic region in 
two-impurity Kondo effect can be observed through the 
change of the STS spectra. In particular, the existence 
of the strong antiferromagnetic correlation between lo- 
calized spins are observed as dip structures in the dl/dV. 
Our results indicate the possibility in STM observation 
with the resolution of magnetic interactions on surface 
system, which would contribute to the realization of 
spintronics. 
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